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Abstract 

We consider the Bose-Einstein condensation (BEC) in a relativistic pion gas. The thermodynamic 
limit when the system volume V goes to infinity as well as the role of finite size effects are studied. 
At y — > oo the scaled variance for particle number fluctuations, uj = {AN'^)/{N), converges to finite 
values in the normal phase above the BEC temperature, T > Tq. It diverges as uj (x V^^^ at the BEC 
line T = Tc, and uj (xV at T < Tq in a phase with the BE condensate. Possible experimental signals 
of the pion BEC in finite systems created in high energy proton-proton collisions are discussed. 
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I. INTRODUCTION 



Pions are spin- zero bosons. They are the hghtest hadrons copiously produced in high energy 
colhsions. There were several suggestions to search for the Bose-Einstein condensation (BEC) 

n 

of TT-mesons (see, e.g., Ref. However, no clear experimental signals were found up to 

now. Most of previous proposals of the pion BEC signals were based on an increase of the 
pion momentum spectra in the low (transverse) momentum region. These signals appear to be 
rather weak and they are contaminated by resonance decays to pions. In our recent paper 
it was suggested that the pion number fluctuations strongly increase and may give a prominent 
signal of approaching the BEC. This can be achieved by selecting special samples of collision 
events with high pion multiplicities. 

In the present paper we study the dependence of different physical quantities on the system 
volume V and, in particular, their behavior a.tV-^ oo. As any other phase transition, the BEC 
phase transition has a mathematical meaning in the thermodynamic limit (TL) V ^ oo. To 
define rigorously this limit one needs to start with a finite volume system. Besides, the finite size 
effects are important for an experimental search of the pion BEC fluctuation effects proposed 
in Ref. The size of the pion number fluctuations in the region of the BEC is restricted 
by a finite system volume V. To be definite and taking in mind the physical applications, we 
consider the ideal pion gas. However, the obtained results are more general and can be also 
applied to other Bose gases. 

The paper is organized as follows. In Section II we consider the BEC in the TL. Here we 
emphasize some specific effects of the BEC in relativistic gases. Section HI presents a systematic 
study of the finite size effects for average quantities and for particle number fluctuations. In 
Section IV we discuss the finite size restrictions on the proposed fluctuation signals of the BEC 
in high energy collisions with large pion multiplicities. A summary, presented in Section V, 
closes the paper. 
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II. BEC IN THERMODYNAMIC LIMIT 



A. Phase Diagram 

We consider the relativistic ideal gas of pions. The occupation numbers, Upj , of single 
quantum states, labelled by 3-momenta p, are equal to Upj = 0, 1, . . . , oo , where index j 
enumerates 3 isospin pion states, 7r+, vr^, and tt^. The grand canonical ensemble (GCE) average 
values, fluctuations, and correlations are the following j^, 

1 



exp 



r 



1 



((^^Pj) ) = ((^P,i - i^pj)) ) = (^pj) (1 + (^pj)) = Vp,i , (Ar2pjAnk,i) = t^pj^pk^ji , (2) 

where the relativistic energy of one-particle states is taken as ep = (p^ + m^)^/^ with 
m = 140 MeV being the pion mass (we neglect a small difference between the masses of charged 
and neutral pions), T is the system temperature, and chemical potentials are /i+ = /i + /ig, 
/i_ = /i — fiQ, and fiQ = fi, for 7r+, vr^, and 7r°, respectively. In Eq. ([1]) there are two chemical 
potentials: fiq regulates an electric charge, and /i a total number of pions. In this paper we 
follow our proposal of Ref . and discuss a pion system with jJ,Q= ^ . This corresponds to zero 
electric charge Q which is defined by a difference between the number of tt^ and vr" mesons, 
Q = — = 0. The total pion number density is equal to: 

Ep,,(^p,.) ^ 3 r p^dp 



p{T, /i) = P+ + p+- po 



V 27r2 Jo 



exp 



(^a/p^ + - p^ /T 



- 1 



STm 1 

= p*(r,/i) = y]- K2{nm/T) exp{np/T) , (3) 

n=l 

where pj = {Nj)/V, with j = +, — , 0, are the pion number densities, V is the system volume, 
and K2 is the modified Hankel function. In the TL, i.e. for V 00, the sum over momentum 
states is transformed into the momentum integral, . . . = {V/27c'^) . . .p^dp . The particle 
number density p depends on T and p, and volume V does not enter in Eq. Q. This is 
only valid ai p < m. The number of particles at zero momentum level is then finite, and its 



^ The BEC in the relativistic gas of 'positive' and 'negative' particles at fiq m has been discussed in 
Refs. HflflBH 
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contribution to particle number density goes to zero in the TL. The inequahty /i < m is a 
general restriction in the relativistic Bose gas, and fi = m corresponds to the BEC. The Eq. ([3]) 
gives the following relation between the BEC temperature Tq and total pion number density p 

2 °o 



P = ^'^^T y^--^2 (nm/Tc) exp{nm/Tc) • (4) 

n=l 

A phase diagram of the ideal pion gas in the p — T plane is presented in Fig. [H 
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FIG. 1: The phase diagram of the relativistic ideal pion gas with zero electric charge density pQ = 
pj^ — = 0. The solid line shows the relation T = Tc @ of the BEC. Above this line, T > Tq and 
there is a normal phase described by Eq. ([3]). Under this line, T < Tc and there is a phase with the 
BE condensate described by Eq. ([8]). 



The line of the BEC phase transition is defined by Eq. (jl]), and it is shown by the solid line in 
Fig. [TJ In the non-relativistic limit, Tc/m<^l, using K2{x) = ^J-n j (2x) exp(— x) at x ^ 1, one 
finds from Eq. (jlj), 

Tc = 2tt [3C(3/2)]-2/=^ m"^ p^/^ = 1.592 m"^ p2/^ (5) 
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whereas in the ultra-relativistic hmit^, Tc/m^l, one uses K2{x) = at x <^ 1, and Eq. (jlj) 
gives Q, 

Tc = [3C(3)/7r2]"^/^ = 1.399 p^^^ . (6) 



In Eqs. PED, Cik) = T.T=i'^~^ is the Riemann zeta function ((3/2) = 2.612 and C(3) = 
1.202. The equation ^ corresponds to the well known non-relativistic result (see, e.g., Ref. j^l) 
with pion mass m and 'degeneracy factor' 3. 

The particle number density is inversely proportional to the proper particle volume, p oc r~^. 
Then it follows from Eq. ([5]) for a ratio of the BEC temperature in the atomic gases, Tc{A), to 
that in the pion gas, Tc(7r), in non-relativistic approximation, 

Tc(7r) 



^10^° 



(7) 



Tc{A) m yrj^ J m 
where ta = 10~^ cm and r^r = 10^^'^ cm are typical radiuses of atom and pion, respectively, 
and rriA is the mass of an atom. The Eq. ([7]) shows that Tc(vr) ^ Tc{A), and this happens due 
to <^ rA- 

The equation ([3]) gives the total pion number density at V — > oo in the normal phase T > Tc 
without BE condensate. At T < Tc the total pion number density becomes a sum of two terms. 



P 



PC + p*{T,p = m) 



(8) 



The second term in the r.h.s. of Eq. ([8]) is given by Eq. ([3]). The BE condensate pc defined 
by Eq. (JHl) corresponds to a macroscopic (proportional to V) number of particles at the lowest 
quantum level p = 0. 

To obtain the asymptotic expansion of p*{T,p) given by Eq. ([3]) at /i ^ m — we use the 
identity [exp(y) — 1]^^ = [cth(?//2) — l]/2 and variable substitution, p = \/2m}/'^{m — pY^'^x, 
similar to those in a non-relativistic gas [l^. Then one finds. 



p*(T,m) - p*(T,/i) 



3m3/2 



(m — pY^'^ / x'^dx 



- cth 

- cth 



-y/2m(m — p)x^ + 



p 



2T 

(m — p){x'^ + 1) 
2T 



'0 

3m3/2 
72^ 



cth 



A/2m(m — p)x'^ + m? — m 



2T 



(m - pf/"^ / x^dx 



cth 



'm — p)x^ 



V2' 



poo 


r 1 1 1 


/ x'^dx 




Jo 


-X^ + 1- 



2T 
V27r 



(9) 



(m - pf/^ . 



^ The BE condensate formed in the uhra-relativistic regime has been considered in Refs. 0, as a dark matter 
candidate in cosmological models. 
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At constant density, p*(T, jj) = p*(T = Tq, p = m), one finds in the TL at T ^ Tc + 0, 



TT 



p*(T = Tc,fJ, = m) 



d p{T, p = m) 



dT 



3Tc m3/2 



{m-pf/^ (10) 



Using Eq. (ITOl) one finds the function p{T) at T ^ Tc + 0, 



m — p(T) 



2n^ \dp{T,p = m) 



9T2 m3 



dT 



T=Tc 



C) 



fill 



In the non-relativistic hmit m/T ^ 1 one finds, p{T,p = m) ^ 3C(3/2) [mT/ (27r)]^/^ similar 
to Eq. ([5]). Using Eq. (fTTj) one then obtains [l^ . 



m - /i(T) ^ 9C^(3/2) 



T-T, 



(12) 



Tc I6n \ Tc 

In the ultra-relativistic hmit, T/m ^ 1, one finds, p{T,p = m) = 3C(3) T^/tt'^, similar to 
Eq. Q). This gives. 



m - p{T) ^ 18 C'(3) fTc 

m 



'-C 



IT" 



T-T, 



c 



c 



(13) 



Thus, p{T) — > m and dp/dT a.t T ^ Tc + 0, both p(T) and dp/dT are continuous 
functions at T = Tc. 



B. Specific Heat at Fixed Particle Number Density 



The standard description of the BEC phase transition in a non-relativistic Bose gas 
is discussed in terms of the specific heat per particle at finite volume, Cy/N (see, e.g.. 



Refs. 



m 
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14|). The relativistic analog of this quantity is: 

cy _ 1 f de\ 
7 " P [dfj^ • 



(14) 



The energy density in the TL equals to: 



e{T,p) = PC m + e*{T,p) = pc m + — 

OT-i2 2 °° 

SI m 



p^dp 



exp 



Pc m ^ 



27r- — ' [ n 

n=l 



'nm 



V T 



m 
2^ 



(nm\ (nm\ 

(-) + /f,(-)j 



a/ + — /ij /T 
/np 



exp 



V T 



(15) 
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where pc = and jj, < m aX T > Tc, while pc > and /i = m at T < Tc- The high and 
low temperature behavior of cy/p can be easily found. At T — oo and fixed p, both the Bose 
effects and particle mass become inessential. The energy density f|T5l) behaves as, e = 3Tp, 
thus, cy/p = 3. Note that in a non-relativistic gas at ^ T ^ m, one finds e = (m + 3T/2)p. 
Thus, the 'high-temperature non-relativistic limit' would give cy/p = 3/2. At Tc ^ T —>■ 
the behavior of e at fixed total particle number density, p = pc + p*iT, p = m), is given by. 
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FIG. 2: The solid lines demonstrate the temperature dependence of cy/p (I14p in the XL at fixed 
values of /> = 0.005 and 0.15 fm^'^. The cy/p in the finite system is described by Eq. ()34p (see next 
Section). The dashed lines correspond to V = 10^ fm^ and the dashed-dotted lines to V = 10^ fm^, 
respectively. 

The Fig. [2] shows the temperature dependence of cy/p at fixed p. As seen from Fig. [2l 
cy/p (fT4l) has a maximum at T = Tc- The cy/p is a continuous function of T, whereas its 
temperature derivative has a discontinuity at T = Tc- This discontinuity emerges in the TL 
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V ^ oo and can be classified as a 3rd order phase transition. To estimate the value of cy/p 
(HM aX T = Tc we start from T > Tc when the contribution from p = level to cy/p ( |T4l) 
equals to zero in the TL, and then consider the limit T ^ Tc + 0. We discuss separately the 

non-relativistic and ultra-relativistic approximations. 

Using the asymptotic, K,y{x) 
a non-relativistic limit of Eqs. ( 

rr\ 3/2 



A/7r/(2x) exp(-2;)[l + [Au"^ - l)/8x] at a; > 1 [llj, one finds 
and (fTSl) . respectively, 



p(T,p) = 3 
£(T,/i) = 3 



2tt 



Li 



3/2 1 



15 T 



m 



Li 



5/2 1 



mT 



3/2 



27T 

^^3/2(2) + ^ — Li^/2{z) 
8m 



pm 



(19) 



where z = exp[— (m — /i)/T] and Lik{z) = ^^=iZ^/n^ is the polylogarithm function [15i] . 
At T = Tc it follows, ;z = 1 and Lik{l) = ({k). Using Eqs. ^EEB, and ([I2D one finds at 
T = Tc < m, 

^ 15 C(5/2) 



1.926 . 



py^^^^ 4 C(3/2) 

Using the asymptotic expansion, Ki,{x) = |r(z/)(x/2)^'' at x ^ 1 
in the ultra-relativistic limit T > Tc ^ m, 

14 00 



(20) 



ll|, one finds from Eq. (ITS!) 



£(T,p) 
From Eqs. (1211) and 



3^4 



n=l 



np\ 
T / 



9C(4) 



Ji4 



9C(3) ^3 



p 



(21) 



it then follows at m ^ T = T( 



Cy_ 
P 



12 



C(4) 



10.805 . 



(22) 



T=Tc ^^'^^ 

Different pion number densities correspond to different values of the BEC temperature Tc- The 
Eqs. ( I20|22l) show that cy/p goes at T ^ 00 to its limiting value 3 from below, if Tc is 'small', 
and from above, if Tc is 'large' (see Fig. [2]). 

Note that at the BEC in atomic gases the number of atoms is conserved. Thus, the temper- 
ature dependence of cy/p for the system of atoms at fixed p can be straightforwardly measured. 
This is much more difficult for the pion gas. There is no conservation law of the number of 
pions, and the special experimental procedure is needed to form the statistical ensemble with 
fixed number of pions. 
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III. FINITE SIZE EFFECTS 



A. Chemical Potential at Finite Volume 

The standard introduction of pc with Eq. ([8]) is rather formal. To have a more reahstic 
picture, one needs to start with finite volume system and consider the limit V ^ oo explicitly. 
The main problem is that the substitution, X]p • • • — (^/^^r^) . . . p^dp, becomes invalid 
below the BEC line. We consider separately the contribution to the total pion density from 
the two lower quantum states, 

1 , 3 1 3 6 



27r 



oo 



+ — / P'dP ^ X • (23) 



pi exp 



a/ + p"^ — jjj /T 



The first term in the r.h.s. of Eq. fl23|) corresponds to the lowest momentum level p = 0, 
the second one to the first excited level pi = 27rV^^^^ with the degeneracy factor 6, and the 
third term approximates the contribution from levels with p > pi = 2txV^^I'^ . Note that this 



corresponds to free particles in a box with periodic boundary conditions (see, e.g., Ref 161]). 
At any finite V the equality fi = m is forbidden as it would lead to the infinite value of particle 
number density at p = level. 

At T < Tc in the TL one expects a finite non-zero particle density pc at the p = level. 
This requires {m — p) /T = 5 (x V^^ at ^ ^ oo. The particle number density at the p = pi 
level can be then estimated as, 

18 V-^ 18 V-^ V-^ 



Pi 



exp 



oc 



V-^'"" , (24) 



and it goes to zero aXV oo. Thus, the second term in the r.h.s. of Eq. fl23|) can be neglected 
in the TL. One can also extend the lower limit of integration in the third term in the r.h.s. of 
Eq. (l23l) to p = 0, as the region [0,pi] contributes as V^^^"^ ^ in the TL and can be safely 
neglected. Therefore, we consider the pion number density and energy density at large but 
finite V in the following form: 

p = ^ ^ ^— + p*(T,/i) , (25) 

V exp [(m — p) /TJ — 1 

3 Tfi 

e - - ^ — + e\T,p) . (26) 

V exp [[m — p) /TJ — 1 
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Thus, at large V, the zero momentum level defines completely the finite size effects of the pion 
system. 

The behavior of p*(T, /i) at fi —>■ m can be found from Eq. (Q. At large V, Eq. fl25|) takes 
then the following form, 

P = + P*(T,p) = ^ + P*(r,/i = m) - -|- {mTr/'V~6. (27) 

The Eq. ( 1271) can be written as, 

A6^/^ + B6 - 1 = , (28) 

where 

A = {mTf^ = a{T) V , (29) 

V 277 

B = \[p - P*{T,^^ = m)] ^ h{T)V . (30) 

The Eq. ( I28ll for 5 has two complex roots and one real root. An asymptotic behavior at K ^ cxd 
of the physical (real) root can be easily found. KtT <Tc it follows from Eq. (l30l) that h{T) > 0, 
and one finds from Eq. (1281) at large V , 

6 ^ Iv-^ (31) 





From Eq. ( l30l) one finds that 6 = at T = Tc. In this case, Eq. ( 1281) gives, 

5 - 4? • (32) 

The Eq. f l28|) can be also used at T > Tc if T is close to T^, thus, 5^1. In this case it follows 
from Eq. fl30|) that 6(T) < 0, and one finds from Eq. fl28l) . 

5 = ^ . (33) 

Thus, 5 is small but finite at — oo, and /i remains smaller than m in the TL. The temperature 
dependence of chemical potential fi = ^i{T) for V = 10^ fm^ and 10^ fm^ at fixed pion number 
density p is shown in Fig. [31 

The value of cy/p ([141) at finite volume V is calculated as, 

cv _ I f de\ 3 exp [(m — /i) /T] m {m — ji + jj,' T) 

J ^ ~p [dfj^ y ^ JV (exp [(m - /x) /T] - 1)' 



3 roo exp 



^A/m^ + P"^ — P^ /T a/^^ + p2 ^^^2 _|_ p2 _ ^ _|_ ^' 2^ 



^""'^ "^0 (^exp [(^v^m2+p2 _ /t] - 
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FIG. 3: The chemical potential /i as a function of temperature T at fixed particle number density p. 
The solid line presents the behavior in the TL V ^ oo. The dashed line corresponds to y = 10^ fm^, 
and dashed-dotted line to y = 10'^ fm'^. The vertical dotted line indicates the BEG temperature Tc- 
The left panel corresponds to p = 0.05 fm~'^, the right one to p = 0.15 fm~^. 

where /i' = {dfi/dT)py. The temperature dependence of cy/p flMl) at several fixed values of V 
is shown in Fig. [2] by the dashed and dashed-dotted lines. 

B. Particle Number Fluctuations 

The variance of particle number fluctuations in the GCE at finite V is: 

3 



{AN') ^ {{N - {N)f) = + 



exp [(m — jj) /T\ — 1 



+ 



p'^dp 



exp 



y^mM-p^ — fi) /T 



{exp [(m - /i) /T] - 1} 27r ^ jexp ^(^^m^ + p2 - /t] - l} 



2 ' 



(35) 



where the first two terms in the r.h.s. of Eq. fl35|) correspond to particles of the lowest level 
p = 0, and the third term to particles with p > 0. We will use the scaled variance, 

{AN') 



UJ 



{N) 



(36) 



as the measure of particle number fluctuations. The numerical results for the scaled variance 
are shown in Fig. HI At T > Tc the parameter 6 goes to the finite limit (l33l) at V ^ oo. This 
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FIG. 4: The dashed Unes show the GCE scaled variance (j36p for the pion gas as a function of 
temperature T for V = 10^ fm'^, 10^ fm^, 10^ fm^ (from top to bottom). The vertical dotted line 
indicates the BEG temperature Tc- The solid line shows the uj (j36|) in the TL V ^ oo. The left panel 
corresponds to p = 0.05 fm~^, the right one to p = 0.15 fm~'^ . 



leads to the finite value of u ( ]36l) in the TL. At T < Tc one finds from Eq. (]35l) in the TL, 

{AN^} ^ 35-^ + - V a d'^''' . (37) 

This gives, 

^ - ^ = 3 p'^ + ' (38) 

where a = a(T) is defined in Eq. ( 12^ . The substitution of 6 in Eq. (155]) from (^11) . gives for 
T < Tc and oo, 

u; ^ 3b^ p-'V + h"^ V^'^ = ujc + UJ* . (39) 

The uoc in the r.h.s. of Eq. ( l39l) is proportional to V and corresponds to the particle number 
fluctuations in the BE condensate, i.e. at the p = level, uc = J2j{i^^p=o,jY) / {^) ■ The 
second term, u* is proportional to V^^"^. It comes from the fluctuation of particle numbers at 
p > levels, UJ* = Epj;p>o((^^p,i)^)/(^)- At T ^ 0, one finds a ^ and 6 ^ p/3. This 
gives the maximal value of the scaled variance, u = pV/3 = {N)/3, for given p and V values. 
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FIG. 5: The upper panel shows the ratio of condensate particle number density to the total particle 
number density, pc/Pi as functions of T for V = 10^, 10^, 10^ fm'^, and in the TL V — > oo. The 
lower panel shows the ratio of particle number fluctuations in condensate to the total particle number 
fluctuations, ujc/uj, as functions of T for the same volumes. The vertical dotted line indicates the 
BEG temperature Tc- The left panel corresponds to p = 0.05 fm^^, the right one to p = 0.15 fm^^ . 

The substitution of 6 in Eq. (155]) from (1521) . gives for T = Tc and V ^ oo, 

u ^ 3 + 3 ^4/3 ^1/3 ^ ^ ^ 

The Fig. [5] demonstrates the ratios pc/p and uc/uj as the functions of T for = 10^, 10^, 
10"^ fm^, and at \^ ^ oo. In the TL V — > oo, one finds pc at T > Tc- The value of 
Pc starts to increase from zero at T = to p at T ^ 0. Thus, pc remains a continuous 
function of T in the TL. In contrast to this, both ujc and uj* have discontinuities at T = Tc. 
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They both go to infinity in the TL V oo. The ujc/uj ratio equals to zero at T > Tq, 'jumps' 
from to 2/3 at T = Tc, and further continuously approaches to 1 at T ^ 0. At T = Tc 
the contribution of p = level to particle density, pc, is negligible at V ^ oo, but the scaled 
variance uc from this level equals 2/3 of the total scaled variance u and diverges as V^^^. We 
conclude this section by stressing that the particle number fiuctuations expressed by the scaled 
variance uj looks as a very promising quantity to search for the BEC in the pion gas. 



IV. BEC FLUCTUATION SIGNALS IN HIGH MULTIPLICITY EVENTS 

In the GCE, the scaled variances for different charge pion states, j = +, — ,0, are equal to 
each other and equal to the scaled variance uj for total number of pions, 

Epj(^P,i) 

There is a qualitative difference in the properties of the mean multiplicity and the scaled vari- 
ance of multiplicity distribution in statistical models. In the case of the mean multiplicity results 
obtained with the GCE, canonical ensemble, and micro-canonical ensemble (MCE) approach 
each other in the TL. One refers here to the thermodynamical equivalence of the statistical en- 
sembles. It was recently found , [lol, [2^ that corresponding results for the scaled variance 
are different in different ensembles, and thus the scaled variance is sensitive to conservation 
laws obeyed by a statistical system. The differences are preserved in the thermodynamic limit. 
Therefore, the pion number densities are the same in different statistical ensembles, but this is 
not the case for the scaled variances of pion fluctuations. The pion number fluctuations in the 
system with flxed electric charge, Q = 0, total pion number, A^, and total energy, E, should be 
treated in the MCE. The volume V is one more MCE parameter. 
The MCE microscopic correlators equal to (see also Refs. 

,Q): 



2 2 



QjQi A(e2) + epek A{n^) - {e^ + ek)A(7re) 



(42) 



_A(g2) A(7r2)A(e2) - (A(7re))2 

where q+ = 1, = -1, go = , A{q^) = Ep,j<lK,j > = Epj^h , ^(e') = 

Epj ^p^pJ ' ^(^^) — X]p j ^p'^^p.i • Note that the flrst term in the r.h.s. of Eq. (H2l) corresponds 
to the GCE ([2]). From Eq. (H2|) one notices that the MCE fluctuations of each mode p are 
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reduced, and the (anti) correlations between different modes p 7^ k and between different charge 
states appear. This results in a suppression of scaled variance Umce in a comparison with the 
corresponding one uj in the GCE. Note that the MCE microscopic correlators fH2l) . although 
being different from that in the GCE, are expressed with the quantities calculated in the GCE. 
The straightforward calculations lead to the following MCE scaled variance for vr'^-mesons ^|: 

Ep,k(^^P,o Ank,o)mce _ 2 

Due to conditions, = iV_ and + + A^o = N, it follows, 00^^^ = uj^^J^ = u/Q and 
<ce = = ^/3, where iV,, = iV+ + N_. 

The pion number fluctuations can be studied in high energy particle and/or nuclei collisions. 
To search for the BEG fluctuation signals one needs the event-by-event identifications of both 
charge and neutral pions. Unfortunately, in most event-by-event studies, only charge pions are 
detected. In this case the global conservation laws would lead to the strong suppression of the 
particle number fluctuations, see also Ref. j^, and no anomalous BEG fluctuations would be 
seen. 

As an example we consider the high vr-multiplicity events in p+p collisions at the beam energy 
of 70 GeV (see Ref. 2l|). In the reaction p + p p + p + N with small final proton momenta in 



the c.m.s., the total cm. energy of created pions is E = 2mp = 9.7 GeV. The estimates 22 1 
reveal a possibility to accumulate the samples of events with fixed = 30 ^ 50 and have the 
full pion identification. Note that for this reaction the kinematic limit is N"^""^ = E/itlt^ = 69. 
To define the MCE pion system one needs to assume the value of V, in addition to given fixed 
values of Q = 0, E = 9.7 GeV, and A^. The T and fi parameters of the GCE can be then 
estimated from the following equations, 

E = Ve{T,f^;V) , A^ = V piT,fx;V) . (44) 

In calculating the e and p in Eq. (jUj) we take into account the finite volume effects according 
to Eqs. fl25ti26p as it is discussed in Sec. II. Several 'trajectories' with fixed energy density are 
shown in Fig. [6] starting from the line /i = in the pion gas in the p — T phase diagram. The 
MCE scaled variance of vr*^ number fluctuations, uj^^^, increases with increasing of A^. The 
maximal value it reaches at T — 0, 

2 2 / T\Tmax\ 

- 3 (1 + = 3 + ^r) - ^^^^ 



15 




P (fin ') 

FIG. 6: The phase diagram of the ideal pion gas with zero net electric charge. The dashed line 
corresponds to p = p*{T, fi = 0) and the solid line to the BEG T = Tc both calculated in the TL 
V — > oo. The dashed-dotted lines present the trajectories in the p—T plane with fixed energy densities, 
e = 6, 20, 60 MeV/fm^ calculated for the finite pion system with total energy E = 9.7 GeV according 
to Eq. (j26p . The dotted lines show the same trajectories calculated in the TL V ^ oo. The total 
numbers of pions N marked along the dashed-dotted lines correspond to 3 points: p = 0, T = Tc, 
and T = for E = 9.7 GeV. 



In Fig. m, ^nice is shown as the function of N. Different possibilities of fixed energy densities 
and fixed particle number densities are considered. One way or another, an increase of N leads 
to a strong increase of the fluctuations of A^^q and Nch numbers due to the BEC effects. 

The large fluctuations of Nq/Nch = f ratio were also suggested (see, e.g., Ref. [23|) as a 
possible signal for the disoriented chiral condensate (DCC). The DCC leads to the distribution 
of / in the form, dW (f) / df = l/{2y/f). The thermal Bose gas corresponds to the /-distribution 
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FIG. 7: The scaled variance of neutral pions in the MCE is presented as the function of the total 
number of pions A^. Three solid lines correspond to different energy densities, e = 6, 20, 60 MeV/fm^ 
(from bottom to top), calculated according to Eq. (j26p . Two dashed-dotted lines correspond to 
different particle number densities, p = 0.05, 0.15 fm~'^ (from bottom to top), calculated according 
to Eq. (I25|). The scaled variance ujI^^^ is given by Eq. (03]), with uj ([36]) and {/^N"^) ([35]). The total 
energy of the pion system is assumed to be fixed, = 9.7 GeV. 

centered at / = 1/2. Therefore, /-distributions from BEC and DCC are very different, and 
this gives a possibility to distinguish between these two phenomena. 



V. SUMMARY 



The idea for searching the pion BEC as an anomalous increase of the pion number fluctu- 
ations was suggested in our previous paper [2I. The fluctuation signals of the BEC have been 
discussed in Ref. 2! in the thermodynamic limit. At V ^ 00, it follows, uj = 00 ai T < Tc- 



This is evidently not the case for the finite systems. At finite V the scaled variance u) of the 
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pion number fluctuation is finite for all possible combinations of the statistical system param- 
eters. The lo demonstrates different dependence on the system volume V in different parts of 
the p — T phase diagram. In the TL — > oo, it follows that uj converges to a finite value at 
T > Tc- It increases as oc V^^^ at the BEC line T = Tc-, and it is proportional to the system 
volume, cij oc at T < Tc- The statistical model description gives no answer on the value of V 
for given E and A^. The system volume remains a free model parameter. Thus, the statistical 
model does not suggest an exact quantitative predictions for the A^-dependence of cj^ce ^^"^ 
cu^^g in the sample of high energy collision events. However, the qualitative prediction looks 
rather clear: with increasing of N the pion system approaches the conditions of the BEC. One 
observes an anomalous increase of the scaled variances of neutral and charged pion number 
fluctuations. The size of this increase is restricted by the finite size of the pion system. In turn, 
a size of the created pion system (maximal possible values of and V) should increase with 
the collision energy. 
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